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Abstract 

The interaction of a magnetic flux vortex with weak external fields is considered 
in the framework of the Abelian Higgs model. The approach is based on the 
calculation of the zero-mode excitation probability in the external field. The 
excitation of the field configuration is found perturbatively. As an example we 
consider the effect of interaction with an external current. The linear in the scalar 
field perturbation is also considered. 



Introduction 

As is known, the Abelian Higgs model in 2+1 dimensions possesses a finite-energy static 
topological nontrivial solution of the vortex type [1]. In the nonrelativistic limit the 
model has the form of the Ginzburg-Landau theory, which phenomenologically describes 
the magnetic flux vortices in the type II superconductor [2]. On the other hand, the 
extrapolation of the vortex solution of the Abelian Higgs model to 3+1 dimensions gives 
U(l) string configuration, which could be produced as a topological defect at the early 
Universe [3] . Very recently a new class of closed vortex ring solutions was discovered in 
3+1 dimensional SU(2) Yang-Mills-Higgs theory [4]. Thus the strings or vortices have 
played many interesting roles in the interplay between high energy physics, condensed 
matter and cosmology. 

The dynamical properties of vortices in the Abelian Higgs model and some of its 
modifications have been considered from different points of view. The much studied 
question seems to be a problem of interaction between strings. 
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The general proof of the existence of multivortex configurations was constructed in 
[5, 6] and a method for obtaining asymptotic multisoliton solutions in gauge theories 
was given in [7]. There is a nice description of vortices dynamics based on the moduli 
space technique [8, 9]. In this scheme low energy soliton dynamics is approximated 
by geodesic motion on the space of the collective coordinates of static multivortices 
configuration with respect to the metric induced by functional of the kinetic energy. 
Since no exact static multivortices solutions are known, some numerical calculations 
have been performed. The numerical investigation of the interaction between well- 
separated vortices was carried out by Shellard and Ruback [10]. An analytical study of 
the interaction between vortices was done recently in [11]. 

The general result of these works is that there are no long-range forces between 
static vortices. There is a difference from BPS multimonopole configuration where the 
repulsive and attractive long-range forces between monopoles exactly compencate, BPS 
monopoles experience no net interaction. That makes the problem of the string dynamics 
in an external field more compicated comparing to the case of monopoles. Actially the 
interaction of vortices with external fields was not much studied. 

There is a regular perturbation scheme used to describe the motion of solitons under 
external force [12]. This approach is based on the calculation of the probability of 
the excitation of corresponding translations zero modes in an external field. In our 
previous publication [12] we discussed the application to the problem of motion of 1 + 1 
dimensional kink in 4 model and to the case of the interaction between the 't Hooft- 
Polyakov monopole and an external weak field. 

In the present note we would like to apply this formalism to the problem of interaction 
between 2+1 dimensional vortex solution in the Abelian Higgs model and an external 
field. 

1 Current-current interaction 

The Lagrangian of the Abelian Higgs model in (2+1) dimensional space-time is given 
(in the units c = 1, H — 1): 

L = \F, V F^ + 1 -{D,4>Y{D^) + ±(\ (f) \z- v *)\ (1) 

where a U(l) real gauge potential = (A ,Ak), k = 1,2, is coupled to a charged 
complex scalar field = 4>i + i<t>2- Here e is the gauge coupling constant, A is the scalar 
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field self-coupling constant, v is the vacuum expectation value of the modulus of the 
scalar field and 

D ll (f> = d^(f>-ieA ll (f>, or D^ a = 8^ a - ee ah A^\ a = 1,2 (2) 

are the field strength tensor and the covariant derivative correspondingly. Hereafter, 
the space indices are m,n,k... = 1,2 and the real components of the complex Higgs 
field are labeled by the indices a,b,c... = 1,2. 

This model is a relativistic analog of the Landau- Ginzburg theory of syperconduc- 
tivity. The finiteness of the energy determines the asymptotic form of the fields. There 
are nontrivial static vortex solutions of this model depending on the polar coordinates 
r and if [1], [2] which are described by the Nielsen-Olesen ansatz: 

77 

A = A r = 0; A v = — (1 - KM) ; = vHMe m,p , where neZ. (3) 

er 

Here the boundary conditions on the structure functions K(r), H(r) are: 

K(r) ™0, H(r) ™ 1; 

K( r ) ^ 1, H{r) ^ 

and the only nonzero component of the field strength tensor is F xy = —{n/er)K'{r) = 
B z . 

Let us consider the interaction of this configuration with an external field. Note that 
because both vector and scalar fields are coupled in the vortex configuration, there are 
a few different ways to introduce such an interaction. For example, one could add to 
the Lagrangian (1) an extra current-current perturbation term that is linear in external 
current J m as well as in current of scalar field 

L\Z = J m J m = - W{D m cf>) - {D m 4>y<P\ J m . (4) 
As a result the classical field equations take the form 

d m F mn = |[0*( J D„0)-( J D n 0)*0]+^ 1) ; 
D m D m< p a = A(0V'-^ 2 )0 a + ^ a(2) (5) 

where the last terms represents the external force acting on the configuration. They 
read 
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J™ = e 2 |0| 2 J„; = 2e e ah D n $J n 



(6) 



The effect of such a perturbation is that there are corrections to the static vortex 
solution (3). In an analogy with the case of two-dimensional A0 4 model [12] these 
corrections could be expanded in powers of external perturbation: 



where (A m ) , (0 a )o correspond to the classical n-vortices solution given by the zeroth- 
order approximation (3) and the fluctuations of the vector and scalar fields on this 
background are of the same order as perturbation J m . 

To the first order corrections, they can be found from the equations describing the 
fluctuations of vector and scalar fields: 



where we have decomposed the complex field via x = Xi + and the background field 
gauge d m a m = ie{<f)*x - X*<P) is used. 

Now one can apply the same approach that was already used in [12], i.e. write the 
expansion of the fields a m (r,t), x a (r, t) on the eigenmodes of the matrix T> 2 of second 
functional derivatives of the action with respect to the fields A m , <p a 



■m 



(A m ) + a m + . . . ; 



r = (0 a ) o + x a + • • • 




(7) 




) 



(d 2 m - eWK 2ee ab D n rx b 
-2ee ab D m (f) b a m D m D m x a 



In matrix notation the equation of motion (7) can be rewritten in the form 




(8) 



where 
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We seek for the solution of Eq. (8) in the form of an expansion 

00 

/(r,0 = E^(00(r) (9) 



i=0 



on the complete set of eigenfunctions d(r) of the operator V 2 . These eigenfunctions 
consist of a vector and a scalar component: d(r) = ^ ^^a^) ^ describing the fluctu- 
ations of the corresponding fields on the vortex background [5], [8]. Thus, there are 
indices of two kind: index i is the number of a mode and the indices m,a correspond to 
its spatial and 'isotopic' components. 

The substitution of expansion (9) into Eq.(7) results in the following system of 
equations for coefficients Cj(i): 

00 

E (pi + tfa) rj mi (r) - 2ee abX b D m r = ^(r); 

ufd) 7j?(r) - 2ee ab a m D m <p b = F°<?\t), (10) 



=0 

00 



E(A + 



i=0 



This is a system describing two sets of coupled forced oscillators. 

It is known that among all fluctuations on the vortex background there are modes 
with zero energy (fl = uj = 0) [5]. Such modes are collective translation coordinates 
of the multivortices configuration. All the other solutions of the system (10) correspond 
to the oscillations on the classical background. Unfortunately neither analitical solution 
for the vortex structure functions nor the eigenfunctions ?7 m j(r), rjf (r) are known. Nev- 
ertheless, some information about the vortex dynamic can be obtained from the system 
(10). 

Following the approach by Manton [13], one can find the acceleration of the vortex 
under perturbation J m if the excitations of the zero modes are treated as a nontrivial 
time dependent translation of the n- vortex configuration (3). 

However, the structure of the Lagrangian (1) suggests [14] that the normalizable zero 
modes are not only translations of the topologically non-trivial configurations but 

C'« WS C„(r)=(^) 

where 

^(r) (fc) = F kn = d k A n - d n A k ; ?f(r) (fc) = D k <f = d k <p a - ee ab <p b A k . (11) 
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Here the index k corresponds to the translation in the direction f k . 

Similarly to the case of 3+1 dimensional Georgi-Glashow model, these normalizable 
zero modes coincide with the pure translational quasi-zero modes of the vector and 
scalar fields rj n {r)^ = d k A n ; rj a (r)^ = d k <p a up to a gauge transformation with a 
special choice of the parameter which is just the gauge potential A k itself. These modes 
are normalized in such a way that makes Cq in expansion (9) equal to the displacement 
of the vortex in the direction f k : 

A n (r + St) « A n (r) + d k A n (r)5x k = A n (r) + C (t)^(r) (fc); 

a (r + Sr) » a (r) + d k( p a (r)5x k = a (r) + C {t)rj a (r)<*> (12) 

Now we can project the Eq.(lO) onto the zero modes (11) which yield the equation: 

C J d 2 x [(r/») 2 + (Vm(r)) 2 } ~ 2e J d 2 x e ah {D^x'Vm + a m D m 4> h V a } 
= J d 2 xF m (r)r) m (r)W + J d 2 xT a (r)r) a (vf\ (13) 

The second term on the left-hand side of this equation describes the transitions 
between the vortex zero modes i] a , r\ m and all the other fluctuations 77", r\ mi on the vortex 
background. A substitution of the expansion (9) gives 

- 2e J d 2 x e ab {D m <f) a x b V m + a m D m ^ h ri a> \ 

= 2e f d 2 xe ab ^C l {t)U {m \\r lm + r 1 a r 1 h{m) r 1 A (14) 

The first term in the sum (14) is equal to zero: the vortex configuration behave as 
a particle-like object and there is no effect of mutual transitions between the collective 
coordinates of scalar and vector fields. All the other terms ii 7^ 0) describe the effect 
of a bremsstrahlung of both vector and scalar massive fields from a vortex accelerated 
by an external force. Such effects of radiation are suppressed as ~ exp{— m s ^r} and, if 
scalar (vector) fields fluctuation are very heavy, they can be considered as an additional 
small perturbation of the second order. But if m s ^ are small, the contribution of this 
transition term (14) could be of the same order as the probability of the zero mode 
excitation. We will not consider here this situation and neglect the contribution from 
the term (14). 

Note, that the first integral on the left-hand side of the Eq.(13) gives a very simple 
result. Because of the virial theorem (see e.g, [15]), the kinetic energy of a vortex is 
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equal to the potential energy and we obtain 



/ d 2 x [(^(r)) 2 + ( Vm (r)) 2 } = J d 2 x (Ft + (D k( f>) 2 ) 
= f<?x{±(F& + (D k <l,)*)+V[<l>]} = M 



(15) 



that is simple the energy of static vortex per unit length, or its mass. Note that the same 
relation between the normalization factors of monopole zero modes and the monopole 
mass holds in 't Hooft-Polyakov model [12]. 

Now we turn to the right-hand side of the Eq.(13). Suppose that the external 
constant current J, m is directed along the x axis: J m = (J, 0). Taking into account the 
definition (6) of external force acting on a vortex and substitute the ansatz (3), one 
could find 



h = f rf 2 ^i 1 )(r)7 ?m (r) (fe) = e 2 J m J dxdy <py a F km 

/ 2 = |rf 2 x^ a(2) (r)^(r) (fe) = -2eJ m j ' dxdy e ab D m <p a . (16) 

Thus the non-trivial result gives the projection onto the zero mode component C^( r ) : 

h = e 2 J J dxdy (f) a <f) a F X y = 27rev 2 Jn J drH 2 {r) dK ^ ■ 



dr 

I 2 = -2eJ J dxdy e ab D x <p a D y( j) b = -Aixev 2 Jn J drH(r)^p- 

= -2irev 2 Jn, (17) 

i.e. the only zero modes orthogonal to the external force are excited 1 

Obviously, the result depends upon the relation between the masses of scalar (m 2 = 
2\v 2 ) and vector (m 2 = e 2 v 2 ) particles. For example, in the London limit m s ^> m v , 
and on the distance ranges at mj 1 C r < m~ x one could neglect the core structure of 
scalar field, i.e. suppose that H ~ 1 everywhere. Then we have for the second integral 
in (16) 

2nev 2 Jn J drH 2 (r) dK ^ ~ 2nev 2 Jn J dr dK ^ = 2i\ev 2 Jn. 



dr 



It is interesting to compare this conclusion with another problem of the interaction between well 
separated vortices discussed in [16], where the intervortex forces lead to the well known effect of 7r/2- 
scattering of two vortices by head-on collision [8] . 
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Thus, taking into account the definition of the magnetic flux <3> n = / d 2 xF xy = — , 
we finally obtain 

MC = -4irev 2 Jn = -2v 2 e 2 J$ n = -2m 2 v J<$> n (18) 
and the acceleration of a vortex along the direction orthogonal to the external force is 

W= 2J ^$n (19) 

2 A linear perturbation of the scalar field 

Alongside with the Lagrangian of interaction (4) there is another possibility to intro- 
duce linear on scalar field interaction between the vortex scalar field and and external 
perturbation: 

L\ nt = e V - (cpe~ imLp + <f)*e imip ) = eve ab <p a n b = ev 2 H{r) cos(n - m)ip, (20) 

where e <C 1 is a perturbation parameter and n a = (— sinmip, cosmip) is a unit vector. 
Obviously this is a perturbation of the vortex configuration connected with the external 
scalar field of another vortex of topological charge m. The case m = corresponds to 
the topologically trivial external constant scalar field coupled with only one component 
of the vortex scalar field <f) 2 = vH(r) cos tp. 

We will see that the effect of such an interaction term depends from the topology 
of external configuration. Indeed, if both vortices have the same magnetic flux (i.e., 
if m — n), the only effect of the additional term (20) is a small increasement of the 
configuration mass. But if, for example n — m = 1, this term lifts degeneration of 
vacuum as it takes place in the case of 2D A0 4 model or 't Hooft-Polyakov monopole 
[12]. 

The difference from the above discussed situation is that the interaction term (4) 
now only affects the scalar component of the vortex configuration. Indeed, the field 
equations are still given by Eq.(13), but the external force acting on the vortex is now 

jrO-) = o ; = £V £abH b (21) 

Projection of this formula onto zero modes along rc-axis gives 

J d 2 xf a {r)ri a {r) {k) = ev J d 2 xe ab D x <f) a n b (22) 
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= ev 2 



ev 2 7i 



J d 2 x jii 7 cos tp cos(n — m)ip + —HK sin ip sin(n — m)<p 

($m,n+l "I - ^m,n— 1 

) / drrH' + (5 m ,n-i 



i.e. the vortices move under an external force described by (20) only if m — n ± 1. On 
the same way one can see that the projection of (21) onto zero modes along y-axis is 
trivial: 

ev J d 2 xe ab D y 4> a n b 

J d 2 x jif ' sin ip cos(n — m)(p HK cos ip sin(n — m)y? j = (23) 



= ev 2 



for any values of m, n. 

In order to estimate the integrals let us consider a limiting case. In the Bogomol'nyi 
limit the first order equations on the shape functions are simply [17] 

H' = -HK; ^K' = m 2 (H 2 -l) 
r r 

Thus the force acting on the vortex in this limit goes to 

T x — > 2irenv 2 5 m ^ 1 J drHK; T y = 0. (24) 

Thus, the effect of interaction is nontrivial only if m — n — 1. This is just the situation 
when the perturbation term, being considered as a correction to the Higgs potential, 
lifts the degeneration of the vacuum 2 . For example we can consider a simple case 
when a vortex of a unitary magnetic flux interacts with a trivial external homogenious 
perturbation. Note, that then we again have the situation when the force is orthogonal 
to the external perturbation. 

Furter aproximation is to neglect the core structure of the vortex, i.e. suppose, as it 
was done at the end of the previous section, that everywhere H ~ 1, K 6(0). Then 
the force acts on the configuration in the London limit and can be approximated as 

2 

sin 

T x — ► nnev 2 = — v -<$> n ] T y = 0, 

e 

2 This corresponds to the so-called thin-wall approximation of the problem of spontaneous vacuum 
decay [18]. 
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that gives the acceleration of the vortex 



This result can be compared with the final formula from the previous section (19). 
In both situations the force is orthogonal to the perturbation and the acceleration is 
proportional to the magnetic flux of the vortex and m%. 



Conclusions 

This work is based on a paper with V. Kiselev [12]. We have found a regular perturbation 
scheme to describe the motion of vortices in the Abelian Higgs model under external 
force. The key point of our approach is to treat the excitation of the zero modes 
of the vortex as a nontrivial time- dependent translation of the whole configuration. 
The amplitude of this excitation can be calculated from the field equation. We have 
considered two different kinds of external perturbation connected with current-current 
interaction and with an external scalar field. In both cases the force acting on the vortex 
is orthogonal to the external perturbation. 

I would like to thank Valera Kiselev for very helpful discussions and many years of 
fruitful collaboration. I am very grateful to D. Maison and V.I. Zacharov for hospitality 
at the Max Planck Institut fur Physik, Werner- Heisenberg Institut, Munich, where this 
work was completed. 
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